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CHAPTER  1 


INTRODUCTION 

An  investigation  of  the  problem  of  a wire  inside  a cavity  which 
is  excited  by  an  external  source  has  been  undertaken.  The  effects  of 
this  external  source  are  coupled  to  the  cavity  interior  and  wire  through 
an  aperture  in  the  cavity  wall.  The  currents  excited  upon  the  wire  and 
the  fields  within  the  cavity  are  to  be  determined.  This  boundary  value 
problem  is  an  idealization  of  a wire  in  some  metal  enclosure.  As 
examples,  the  wire  may  be  inside  the  shielding  or  housing  of  an  elec- 
tronic or  mechanical  unit  on  an  aircraft,  or  it  might  simply  pass  from 
one  metal  partition  to  another  through  a region  which  is  essentially 
empty. 

This  paper  deals  primarily  with  the  problem  formulation  (Chapter 
2)  and  the  consideration  of  the  many  numerical  difficulties  which  are 
encountered  in  obtaining  a solution  (Chapter  4).  In  addition,  an  analogous 
two-dimensional  problem  is  considered  (Chapter  3)  and  some  sample 
numerical  results  are  given  (Chapter  5).  In  Chapter  6,  in  addition  to  a 
summary  of  this  work,  several  possible  extensions  are  discussed.  The 
remainder  of  this  introductory  chapter  is  devoted  to  the  modeling  of  the 
problem  and  its  historical  background. 

Modeling  of  the  Physical  Problem 

In  order  to  model  the  system  of  Interest,  a rectangular  box 
having  an  aperture  in  a side  wall  and  enclosing  a wire  is  considered. 
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The  appropriate  geometry  is  shown  in  Figure  1.  The  problem  is  formu- 
lated in  the  frequency  domain.  It  should  be  noted,  however,  that  given 
the  frequency  domain  solution,  desired  time  domain  quantities  can  be 
obtained  by  numerical  inverse  Fourier  transform. 

As  is  usiial  in  the  investigation  of  complex  problems,  simplifi- 
cations must  be  invoked  to  render  the  problem  tractable.  The  assump- 
tions and  conditions  of  the  cavity/wire  problem  are  summarized  as 
follows: 

1.  The  cavity  is  the  interior  region  of  a perfectly  conducting 
rectangular  box. 

2.  The  material  in  the  box  is  uniform,  linear  and  isotropic. 

3.  The  cavity  is  excited  through  a small  aperture  in  a cavity 
wall  such  that  aperture  dipole  approximations  may  be  used. 

4.  The  wire  is  straight,  circular  and  perfectly  conducting,  and 
is  oriented  perpendicular  to  a side  wall  of  the  cavity. 

5.  The  wire  ends  may  or  may  not  be  in  electrical  contact  with 
the  cavity  walls. 

6.  The  wire  is  thin  at  the  frequency  of  operation  and  thin-wire 
assumptions  can  be  utilized. 

Historical  Background 

In  recent  years,  a great  deal  of  work  has  been  done  on  the 
problem  of  scattering  from  wires  in  free  space  and  efficient  techniques 
have  been  developed  to  handle  them  (Butler  and  Wilton  1975,  Wilton  and 
Butler  1976).  More  recently  problems  for  which  a wire  couples  to  an 
aperture  in  an  infinite  planar  screen  have  been  considered  (Lin,  Curtis 
and  Vincent  1976;  Butler  and  Umashankar  1976;  Seidel  and  Butler  1976). 

Historically,  the  work  done  on  the  rectangular  cavity  problem 
has  been  primarily  concerned  with  formulating  the  dyadic  Green's 
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functions  for  the  potentials  and  fields  for  a rectangular  cavity.  An 
expansion  for  the  electric  field  dyad  was  found  by  Weyl  (1913,  1915). 

It  was  pointed  out  by  Teichmann  (1952)  that  the  expansion  was  incom- 
plete. This  was  later  corrected  by  the  addition  of  another  term 
(Teichmann  and  Wigner  1953).  Several  other  authors  have  also  shown  that 
that  the  additional  term  is  necessary  for  completeness  (Kurokawa  1958; 
Collin  1973;  Tai  1973;  Howard  1974). 

The  dyad  for  the  magnetic  vector  potential  for  the  rectangular 
cavity  is  formulated  in  Morse  and  Feshbach  (1953,  pp.  1849-1851). 

More  recently,  this  dyad  and  the  dyad  for  the  electric  field  were 
formulated  by  Tai  and  Rozenfeld  (1976)  using  the  vector  wave  functions 
L,  M and  N (Hansen  1935)  and  by  Rahmat-Samii  (1975)  who  also  derives 
the  dyad  for  the  electric  vector  potential.  Both  of  these  authors  took 
care  to  insure  completeness  of  the  expansions. 

Recently,  Cheng  and  Chen  (1976)  formulated  the  solution  for  the 
problem  of  a rectangular  cavity-backed  aperture  in  an  infinite  screen. 
However,  at  the  present  time,  no  numerical  data  has  been  presented 
using  this  formulation. 

It  should  be  noted  that  at  this  time  there  appears  to  have  been 
no  work  performed  on  the  cavity  containing  an  interior  scatterer.  This 
is  probably  due  to  the  present  uncertainty  of  computational  methods 
for  evaluating  the  dyadic  Green's  fiinctions  in  the  cavity,  especially 
in  or  near  source  regions. 
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CHAPTER  2 


FORMULATION 

For  the  piirposes  of  this  problem,  consider  a perfectly  conduct- 
ing rectangular  cavity.  One  comer  of  this  cavity  is  located  at  the 
origin  of  a cartesian  coordinate  system  (Figtore  1).  The  dimensions  of 
the  cavity  are  denoted  by  a,  b and  c in  the  x,  y and  z directions, 
respectively.  Within  this  cavity,  there  is  a perfectly  conducting, 
round  thin  wire  of  radius  r(r<<X)  which  is  parallel  to  the  z-axis. 

This  wire  may  or  may  not  be  attached  to  either  or  both  walls  of  the 
cavity . 

One  of  the  walls  of  the  cavity  is  perforated  by  a small  aper- 
ture whose  center  is  located  at  r * (x  , y , z ) . The  exterior  region 
to  which  the  aperture  couples  the  cavity  interior  may  be  of  two  differ- 
ent types.  The  cavity  may  be  located  behind  an  infinite,  perfectly 
conducting,  planar  screen  such  that  the  cavity  wall  containing  the 
aperture  is  a portion  of  the  infinite  screen.  Alternatively,  the 
cavity  may  be  situated  in  a free  space  environment.  In  either  case, 
the  excitation  for  the  problem  is  provided  by  sources  in  the  exterior 
region. 

Finally,  it  is  assumed  that  the  medium  in  both  t.he  interior 
and  exterior  regions  is  homogeneous  and  isotropic  and  is  characterized 
by  C^.y)  where  e can  be  complex  for  a lossy  medium.  It  is  assiimed 
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that  the  problem  is  time  harmonic  with  angular  frequency  o)  and  the 
factor  e^^^  has  been  suppressed  throughout. 

Dyadic  Green’s  Functions 

In  order  to  formulate  an  integral  equation  for  this  problem, 
it  is  necessary  to  know  the  Green's  functions  for  the  potentials  and 
the  fields  within  the  interior,  or  cavity,  region.  These  Green's 
fxmctions  are  dyadic  in  nature  and,  as  one  would  expect,  are  singular 
in  the  source  region. 

One  can  define  the  dyadic  Green's  function  for  the  magnetic 
vector  potential  by 


(V^  + k^)  (2.1a) 

n X (k*!  -f  77)  • » 0 on  S (2.1b) 

where  k is  the  wavenumber  of  the  homogeneous,  isotropic  medium  of  the 
cavity  interior,  T is  the  identity  dyad,  n is  an  inward-directed  unit 
normal  vector  on  S where  S is  the  surface  of  the  cavity.  This  Green's 
dyad  has  been  derived  by  Tai  and  Rozenfeld  (1976)  in  terms  of  the 
vector  wave  functions  L,  M and  N and  is  shown  as  a matrix  in  Table  1. 

Once  5^  has  been  determined,  the  Green's  dyads  for  the  electric 
and  magnetic  fields  due  to  an  electric  current  source  can  be  found. 

They  are  defined  by 

gg  - (k*T  + 77)  • (2.2) 
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for  the  electric  field  and 

» V X (2.3) 

for  the  magnetic  field.  The  matrix  forms  of  and  Sj^  are  given  in 
Table  2.  They  were  obtained  by  simply  operating  upon  as  prescribed 
by  (2.2)  and  (2.3). 

It  should  be  noted  that  this  result  for  agrees  with  that 
derived  by  Tai  and  Rozenfeld  (1976)  directly  using  the  vector  wave 
functions.  It  also  would  agree  with  a similar  result  obtained  by 
Rahmat-Samii  (1975)  if  an  obvious  sign  error  in  that  paper  were 
corrected.  Note  that  this  agreement  is  a most  important  point.  In 
these  two  papers,  the  authors  have  taken  great  care  to  insure  complete- 
ness of  the  expansion  functions  for  S in  the  source  region  of  the 
cavity.  The  agreement  between  their  results  and  that  given  here 
demonstrates  the  fact  that  completeness  is  insured  when  the  problem 
is  formulated  through  use  of  the  potentials  and  the  fields  are  then 
derived  from  those  potentials  . 

It  now  remains  to  determine  the  dyadic  Green's  function  for 
the  electric  vector  potential  and  its  related  field  dyads.  Consider 
the  Green's  dyad  for  the  electric  vector  potential  defined  by 

(7*  +k^)  gp(r,r')  » - ! 6(r-f')  (2.4a) 

n • gp  * 0 

, on  S.  (2.4b) 

n X V X gp  ■ 0 
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Table  2.  Dyadic  Green's  Functions  for  the  Electric  and  Magnetic  Fields  Due  to  an  Electric 


cosk  X sink  x'  sink  y sink  y'  sink  z cosk  z'  , etc. 

V V \r*  \f*  •>  *9  ’ 


Rahmat-Samil  (1975)  has  obtained  a solution  for  gp.  It  is  given  in 
matrix  form  in  Table  1.  Again  the  sign  error  has  been  corrected. 

Now  that  gp  is  determined,  the  dyads  for  the  electric  and 
magnetic  fields  due  to  a magnetic  current  source  can  be  defined  by 

8e  “ ■ ^ Sp 

and 

gj^  - (k*!  ♦ W)  • gp  . 

If  g is  written  in  matrix  form,  it  is  found  that 
gg(r.r')  * - Sj^(r'  ,r) 

where  the  tilde  denotes  the  transpose  of  the  dyad. 

As  a matter  of  notation,  an  upper  case  G denotes  a dyad  due 
to  an  electric  current  source;  similarly  a lower  case  g denotes  a 
dyad  due  to  a magnetic  current  source.  The  subscript  A,  F,  e or  h 
denotes  the  particular  potential  or  field  which  is  given  by  the  dyad. 

Before  proceeding  with  the  formulation  of  the  integral 
equation,  it  is  worthwhile  to  consider  a few  of  the  general  properties 
of  these  dyadic  Green's  functions.  Probably  the  most  apparent  property 
is  that  each  component  of  each  dyad  is  in  itself  a triply  infinite 
Fourier  sum.  It  can  be  seen,  however,  that  any  one  of  the  sums  can  be 
performed  analytically  using  one  of  the  following  relationships: 


(2.6) 


(2.7) 
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Jl  ” 2asinTi'5a  sinhax<sinha(a-x^)  [2. 8a) 


y cosk  X cosk  x' 

m-O  k^+a*  ^ ^ 

X 


SHhH?  coshax^cosha(a-x^)  (2 . 8b) 


oo  k 

L kW  ’ 2"sinh6:a  coshax<sinha(a-x^)  (2. 8c) 

X*® 

1 


00 

1 

m“l 


k^+a^ 

X 


sink  x.cosk  x. 

2 ^ X ^ 


2 sinhaa  sinhax^cosha(a-x>)  (2.3d) 


where 

^ , x^  » minCx.x')  , x^^  » maxCx.x’)  , o ^ x , x'  £ a . 

Equations  (2.8)  are  easily  derived  by  considering  the  ordinary  differ- 
ential equation 


1^*  k^l  g(x,x')  » 5(x-x') 

on  the  interval  (0,a)  with  various  combinations  of  unmixed  Dirichlet 
and  Meumann  boundary  conditions.  The  function  g is  then  obtained  by 
a closed  form  construction  to  obtain  the  right-hand  side  of  (2.8)  and 
also  by  a spectral  expansion  [to  produce  the  left-hand  side  of  (2.3)]  . 
.^n  alternative  technique  to  derive  these  equations  has  been  given  by 
Collin  (1960,  p.  581). 
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Another  important  property  of  these  triple  sums  can  be  seen  if 

1 X 

one  remembers  that  both  sinh  x and  cosh  x approach  ^ j e for  large  x. 
Thus,  if  |r-r' | 0,  any  one  of  the  triple  sums  can  be  reduced  to  a 

double  sum  which  is  exponentially  convergent.  This  demonstrates  that 
outside  the  source  region,  all  components  of  all  the  dyads  converge, 
and  in  fact,  converge  exponentially.  Therefore,  all  dyad  components 
are  uniformly  convergent  outside  the  source  region  (Titchmarsh  1939, 
p.  4). 

Since  it  is  valid  to  differentiate  a series  term  by  term  pro- 
vided the  resulting  series  is  uniformly  convergent  (Titchmarsh  1939, 
p.  37),  the  method  used  to  construct  the  Green’s  dyads  for  the  fields 
using  (2.2),  (2.3),  (2.5)  and  (2.6)  is  valid  outside  the  source  region. 
This  leads  to  the  next  observation  regarding  relative  convergence  of 
the  sums.  Note  that  the  effect  of  a differential  operator  on  each  term 
of  any  one  of  the  sums  is  to  introduce  a multiplicative  factor  of  m,  n, 
or  I in  the  numerator.  This  will  slow  the  rate  of  convergence  of  the 
series.  Thus,  for  [f-r' | f 0,  components  of  5^  and  gp  will  exhibit  the 
most  rapid  convergence,  whereas  5^  and  gj^,  which  are  constructed  using 
the  second  order  differential  operator  W*,  will  exhibit  the  slowest 
convergence . 


Integral  Equation  Formulation 

To  formulate  the  problem,  one  first  uses  the  theories  advanced 
for  small  apertures  by  Bethe  (1944).  This  theory  allows  a small  aper- 
ture, whose  center  is  at  the  point  f^,  in  a perfectly  conducting  screen 
to  be  replaced  by  equivalent  electric  and  magnetic  dipole  moments 
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located  at  r . Figure  2 (a  and  b)  depicts  this  equivalence.  The 
electric  dipole  moment  is  normal  to  the  cavity  wall,  and  the  mag- 
netic dipole  moment  is  tangential  to  the  cavity  wall.  Note  that 
now  the  aperture  has  been  shorted,  and  thus  one  has  in  this  equivalent 
problem  a wire  scatterer  in  a rectangular  cavity  excited  by  the  dipole 
sources  P^  and  P^^^.  (A  more  precise  definition  of  these  dipole  moments 
will  be  considered  later  in  this  chapter.) 

It  is  useful  to  make  use  of  yet  another  equivalent  problem.  By 
using  the  equivalence  principle  (Harrington  1961)  one  C2m  replace  the 
wire  scatterer  in  the  cavity  by  unknown  surface  currents  which  are 
located  upon  a mathematical  cylindrical  surface  which  coincides  with 
the  siiTface  of  the  wire  in  the  original  problem.  This  equivalence  is 
shown  in  Figure  2 (b  and  c) . .Now  one  must  force  the  boundary  condition 
that  the  tangential  electric  field  must  vanish  on  the  cylindrical  sur- 
face. When  this  is  accomplished  the  two  problems  are  equivalent  and 
the  surface  currents  in  the  equivalent  problem  will  be  equal  to  the 
surface  currents  induced  upon  the  wire  in  the  original  problem. 

Note  that  in  this  second  equivalent  problem  one  has  a cavity 
whose  interior  is  entirely  homogeneous  and  isotropic,  and  is  driven 
by  the  unknown  sources  J , P and  P . Thus  the  fields  in  the  cavity 
may  be  obtained  by  simply  taking  the  scalar  products  of  the  appropriate 
dyadic  Green's  functions  and  these  sources  and  integrating  over  the 
volume  of  the  cavity. 

The  total  electric  field  at  the  point  f in  the  cavity  can  be 
broken  into  two  parts  in  the  following  manner: 
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eC?)  = + E^(r) 


(2.9) 


where  is  the  field  produced  by  the  aperture  dipole  moments  and 
and  E^  is  the  field  produced  by  the  surface  current  density 
Thus 

E^(?)  » ^ ‘ 

and 

• r f 

E®(?)  - - P (k^T  + W).  I • jgCr')ds'  (2.11) 

wire 

where  n = Ai/e  is  the  characteristic  impedance  of  the  medium  interior 
to  the  cavity. 

Since  the  wire  scatterer  is  assumed  to  be  thin  (r<<X) , one  can 
use  the  traditional  thin  wire  approximations  to  simplify  C2.ll).  These 
approximations  have  been  used  extensively  and  are  summarized  as  follows; 

1.  The  circumferential  component  of  current  on  the  wire  is 
negligible  and  may  be  assumed  to  vanish. 

2.  The  remaining  axial  component  of  current  has  no  circumferen- 
tial variation. 

3.  There  is  no  current  on  the  end  caps  of  the  wire  and  the  axial 
current  must  go  to  zero  at  the  ends  of  the  wire. 

4.  It  is  sufficient  to  require  only  the  axial  component  of  the 
electric  field  to  vanish  on  the  wire  surface. 

If  these  approximations  are  incorporated  into  (2.11),  and  the 
axial  component  of  (2.9)  is  forced  to  vanish  on  the  wire  surface,  one 
obtains  the  following  integral  equation: 
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(k^T  + W)* 


j 

wire 


5^(f,r')  • wJ^(r')ds' 


(2.12) 


where  r is  on  the  wire  surface  and  w is  a unit  vector  parallel  to  the 
axis  of  the  wire. 

It  is  seen  that  (2.12)  is  an  integro-differential  equation  in 
which  all  three  components  of  the  dyadic  kernel  and  all  nine  components 
of  the  dyadic  operator  W will,  in  general,  couple.  As  will  be  seen 
in  Chapter  4,  the  Green’s  functions  are  extremely  difficult  to  calcu- 
late numerically  in  or  near  the  source  region.  Because  of  this  problem, 
the  feasibility  of  solving  (2.12)  for  the  most  general  case  of  arbi- 
trary wire  orientation  is  questionable.  However,  if  one  considers  the 
case  where  w is  equal  to  one  of  the  three  cartesian  unit  vectors  (let 
that  unit  vector  be  z by  convention)  then  (2.12)  reduces  to  a scalar 
integro-differential  equation  as  follows: 


‘iJ.J 

wire 


(?,?’)  J2(z')cls’ 


+ z 


E^(?)  = 0 
(2.13) 


for  r on  the  wire  surface. 

Note  that  if  (2.13)  can  be  inverted,  a solution  for  J.  will  be 

obtained.  However  it  should  be  remembered  that  contains  P and  P , 

e m 

which  as  of  yet  are  unknown.  Thus  it  remains  to  obtain  additional 
constraints  which  will  uniquely  specify  the  values  of  these  dipole 
moments . 
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Aperture  Equivalent  Dipole  Moments 
Before  attempting  to  determine  the  dipole  moments  and  P^, 
it  is  necessary  to  review  the  basis  of  small  aperture  theory.  This 
summary  is  patterned  after  that  in  a recent  paper  by  Butler  (1976)  and 
the  work  of  Collin  (1960).  Consider  an  infinite  perfectly  conducting 
screen  at  z = 0 which  separates  two  half  spaces  of  the  same  properties 
(y,e).  This  screen  is  perforated  by  a small  aperture  centered  about 
the  point  (0,0,0).  If  the  aperture  is  sufficiently  small  and  f is 
sufficiently  far  from  the  aperture,  then  the  fields  at  r due  to  the 
aperture  can  be  approximated  by  the  radiation  from  an  electric  dipole 
with  moment  P^  and  a magnetic  dipole  with  moment  P^^  located  at  (0,0,0) 
which  radiate  in  the  presence  of  the  unperforated  screen.  This  equiva- 
lence is  illustrated  in  Figure  3. 

The  moments  of  the  electric  and  magnetic  dipoles  for  the  right 
half-space  (z>0),  which  are  located  at  (0,0,0+)  are  given  by 

Pg  . ca^(Ej‘='(Q)  - E^‘'''(0))2  (2.14a) 

and 

^m  * ■ “m  * (2.14b) 

where  (E^*^",  H^^”)  are  the  short  circuit  fields  in  the  left  half-space, 
that  is,  the  fields  in  left  half-space  in  the  presence  of  the  unperfor- 
ated screen.  Similarly,  are  the  short  circuit  fields  in 
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Figure  3.  Small  Aperture/Dipole  Moment  Equivalent  Problems. 


the  right  half-space.  The  electric  polarizability  and  the  magnetic 
polarizability  relate  the  specific  excitation  to  the  moments  for  a 
given  aperture.  Polarizabilities  are  available  in  the  literature 
(Bouwkamp  1954;  DeMeulenaere  and  Van  Bladel  1977;  Collin  1960,  pp.  285- 
302;  Cohn  1951,  1952)  for  several  small  apertures.  Specifically,  for 
an  elliptical  aperture  defined  by  + C/^/w*)  ■ 1 for  2.  > w. 


a * 

e 


a 

m 


XX 


WU  ’ 


1 

3 - E(C)  ’ 


C2.15a) 

(2.15b) 


1 ttA’c 

(J)  ECO  - K(0 


(2. 15c) 


and  all  other  components  of  are  zero.  The  square  of  the  eccentricity 
(C)  is  defined  by 


w ^ 

C • 1 - (x) 


and  K and  E are  the  complete  elliptic  integrals  of  the  first  and  second 
kind,  respectively,  as  defined  by  Abranowitz  and  Stegun  (1965,  p.  Sf'O) . 

It  should  be  noted  that,  where  the  small  aperture  theory  is 
based  upon  an  aperture  coupling  two  half-spaces,  in  the  actual  problem 
of  interest  the  interior  region  is  a rectangular  cavity  and  the  exte- 
rior region  may  or  may  not  be  a half-space. 
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First,  consider  the  exterior  region.  Suppose  the  cavity  is 
behind  an  infinite  screen,  such  that  the  exterior  region  is  actually  a 
half-space.  Then  the  short  circuit  exterior  fields  are  easily  deter- 
mined from  a knowledge  of  the  incident  field  by  application  of  physical 
optics. 

However,  if  the  cavity  is  not  behind  a perfectly  conducting 
screen,  the  problem  becomes  more  difficult.  It  now  becomes  necessary 
to  determine  the  short  circuit  fields  on  the  exterior  surface  of  a 
rectangular  box  scatterer.  This  problem  has  been  solved  numerically 
by  Tsai,  Dudley  and  Wilton  (1974).  Since  the  short  circuit  fields  are 
related  to  the  surface  current  and  charge  by 

T ^ ^SC“ 

J “ n X H 
s 

and 

/V  =sc- 
» e n • E , 

these  values  could  also  be  provided  by  experimental  measurements  of 

•SC*  “SC* 

surface  charge  and  current  densities.  Note  that  (E  , H ) and 
(E  , H ) have  been  defined  for  the  problem  of  interest  to  be  the 
short  circuit  fields  in  the  exterior  and  interior  regions,  respective- 
ly.  For  the  remainder  of  this  paper  it  will  be  assumed  that  (E  , 
are  known. 

Now  consider  the  interior  region  of  the  problem  as  illustrated 
in  Figure  2c.  It  is  readily  seen  that  (£  , H ) will  be  driven  by 

the  surface  currents  on  the  wire.  However,  this  region  is  a cavity, 
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and  the  use  of  yet  another  equivalent  problem  is  necessary  in  order  to 
account  for  this  problem  properly.  If  the  method  of  images  (Collin 
1960)  is  applied  to  the  equivalent  problem  depicted  by  Figure  2c,  a new 
equivalent  problem  can  be  obtained  in  which  the  aperture  dipoles  and 
the  wire  currents  have  been  imaged  in  such  a way  that  there  is  a three- 
dimensional  infinite  array  of  image  sources  in  a half-space.  A two- 
dimensional  cross-section  of  that  array  is  given  in  Figure  4.  This 
means  that  are  produced  by  all  of  the  sources  in  that  half- 

space except  for  the  original  aperture  dipoles  and  located  at  r^^. 

Now  if  it  is  realized  that  these  arrays  of  image  sources  in  the 
half-space  are  equivalent  to  the  original  sources  in  the  cavity  (the 
problem  of  Figure  2c),  one  obtains  the  following  relations  for  the  short 
circuit  fields  in  the  cavity  region: 

“ °e  * Pm 

nn  n 


and 


wire 


Gg  (?,?')J^(f-)ds' 


(2.16a) 


a ^ « 


P 

m 


I 

wire 


z J,(r')ds' 

to 


(2.16b) 


where  n is  the  unit  vector  in  the  wall  of  the  aperture  and  the  symbol 
(”')  over  the  dyads  indicates  that  the  original  or  self  term  of  the 
image  series  (Green's  function  for  half-space)  has  been  deleted.  Thus 
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Figure  4.  Equivalent  Interior  Problem  Which  Accounts  for  Effects  of 
Wall  Reflections  and  Wire  Currents  on  Aperture  Dipoles. 
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and 


- I 


j-jkR 

TtFr“ 


.-jkR 

"Hr" 


(2.17a) 


(2.17b) 


where  R » |f-f  |.  Now  the  desired  dyads  may  be  obtained  using  the 

a 

following  relations: 


2^  - (k*!  * TV)  • 2^  . 


(2.18a) 


^ ^ \ * 


(2.18b) 


gg  * - V X gp  . 


(2.18c) 


and 

gj^  - (k*I  + W)  • gp  . (2.18d) 

The  functions  defined  by  (2.17)  and  (2.18)  will  be  referred  to  as 
deleted  Green's  functions  in  an  effort  to  indicate  the  deletion  of 
the  singular,  free  space  Green's  function. 

One  now  defines  a^  and  ^2  cartesian  unit  vectors 

which  are  tangential  to  the  wall  of  the  aperture  in  such  a fashion 
that  a^  X a^  ■ n.  For  example,  if  the  aperture  were  in  the  y * 0 wall 
of  the  cavity,  then  a^  would  be  z,  a-,  would  be  x and  n would  be  y. 
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Note  that  since  has  only  a normal  component  and  has  only  a 

I tangential  component,  there  are  only  three  non-vanishing  unknovms; 

! P . P and  P . If  new  unknowns  E_.  and  H_  are  introduced  such  that 

e m,  m-  T T 

n 1 2 n 


and 


a,  . 


(2.19b) 


then  by  (2.14) 


Q( 

e 


and 


(2.20a) 


(2.20b) 


(2.20c) 


Now,  if  (2.20)  is  substituted  into  (2.16),  which  in  turn  is  substituted 
into  (2.19)  one  obtains  the  following  equations; 
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nn  n 11  nl  1 22  n2  2 


f he  • 

' . nz 
wire 


(2.21a) 


I G^  (r^,r')  Jj(?')ds'  - H®‘''(f^)  , (2.21b) 


wire 


^“eS  <^^a'V^T  ■ “m  ^h*  ^^a'^a^^T  " ' “m  K 

n e a a T^  ^21  e ^ ^22  a a u 


j I Ch,/vf')-'afr')'>»'  ■ "r''^>  • 


wire 


(2.21c) 


It  should  be  noted  that  the  unknowns  E.j.  , Rj.  and  H^.  are 

n i *2 

related  to  the  unknown  dipole  moments  by  the  aperture  polarizabilities. 
Thus  if  (2.20)  is  substituted  into  (2.10)  one  obtains 


2 • E^(r)  • a G (r.r JE-.  - jkm  g (r.f  )H_ 
® ®zn  ^ ‘1  "'ll  ®zl  ‘1 


jkna  g (r.r  )H_ 
‘"22  z2  ^ ^2 


(2.22) 


This  means  that  (2.13)  contains  the  unknowns  E.j.  , H.^  and  H.^.  in 

n I 2 

addition  to  J . The  three  additional  constraints  provided  by  (2.21) 
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when  solved  simultaneously  with  (2.13)  will  provide  a unique  solution 

for  the  currents  on  the  wire  and  the  quantities  E_  , H_  and  H_  . 

‘n  U ^2 

It  is  helpful  at  this  point  to  attempt  to  describe  physically 
the  various  terms  in  (2.21).  First,  it  is  noted  that  the  terms  con- 
taining the  deleted  Green's  fimctions  account  for  the  fact  that  the 
aperture  dipoles  will  be  affected  by  the  fields  reflected  back  from  the 
cavity  walls.  This  is  apparent  if  it  is  remembered  that  these  terms 
are  the  fields  in  the  aperture  due  to  the  array  of  the  images  of  the 
aperture  dipole.  These  virtual  sources  account  for  all  of  the  reflec- 
tions and  multiple  reflections  from  the  cavity  walls. 

It  can  also  be  noted  that  the  integral  term  in  (2.21)  repre- 
sents the  field  scattered  back  into  the  aperture  by  the  wire  and  thus 
accounts  for  the  coupling  between  the  currents  on  the  wire  and  the 
aperture  dipoles.  If  all  of  these  coupling  effects  are  assumed  to  be 

negligible,  (2.21)  reduces  to  (E.j.  , flj.)  » • If  on  the 

n 

other  hand  these  effects  are  not  neglected,  the  difference  between 

(E.J.  , ilj,)  and  (E^^~»  will  in  some  way  reflect  the  degree  of  to 

n 

which  the  cavity  walls  and  the  wire  scatterer  couple  to  the  aperture. 
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CHAPTER  3 

A PRELIMINARY  TWO-DIMENSIONAL  PROBLEM 

In  an  effort  to  gain  insight  for  the  three-dimensional  cavity 
problem,  it  is  helpful  to  consider  an  approximately  analogous  two- 
dimensional  problem.  For  this  case,  the  kernel  functions  would  also 
be  singular  Fourier  series,  although  of  lower  dimensionality.  Thus 
where  in  the  original  problem  the  sums  were  doubly  and  triply- infinite, 
they  are  one-  and  two-dimensional  in  this  two-dimensional  analog.  In 
addition  to  any  insights  which  this  effort  might  provide  toward  the 
solution  of  the  three-dimensional  cavity,  the  solution  to  this  analo- 
gous problem  may  provide  data  pertinent  to  the  physical  problem  being 
modeled.  For  example,  in  this  simpler  model,  it  is  feasible  to  account 
for  the  actual  coupling  through  a large  aperture  rather  than  use  the 
aperture  dipoles  for  small  apertures. 

Formulation 

The  geometry  of  the  analogous  problem  is  shown  in  Figure  5.  It 
consists  of  an  incident  plane  wave  impinging  upon  a perfectly  conduct- 
ing infinite  planar  screen  perforated  by  an  infinite  z-directed  slot  of 
width  d.  The  slot  is  backed  by  a rectangular  cylindrical  cavity  of 
depth  a and  width  b.  Within  the  cavity  are  L z-directed  thin  wires  of 
various  radii.  The  incident  plane  wave  is  polarized  such  that  E is 
parallel  to  the  slot.  Thus  only  the  H and  H field  components  are 
excited. 
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One  now  formulates  the  problem,  using  the  thin  wire  approxi- 


mations described  previously  in  Chapter  2,  by  expressing  E and  H in 

z y 

terms  of  the  unknown  aperture  electric  field  and  the  unknown  surface 

currents  on  the  wires.  Then  by  forcing  the  boundary  conditions  that 

E and  H be  continuous  through  the  aperture  and  E vanish  on  the  wire 
1 y ‘ 

surfaces,  two  coupled  integral  equations  are  obtained.  They  are 

Ej  (y')  Oc|y-y'|)  * jS^(y,y')]dy' 

S2Cx',y,y')dc' 

*■  £.th  wire 

» - 2 (o,y)  , h<y<h+d  (3.1a) 


> k^) 


L 

I 


and 


h+d 


E“  (y')  S2(x,y,y')dy' 


- jkn  I Jo 

i-i 


S2(x,x' ,y,y')dc' 


ith  wire 


0 , (-x,y)  on  surface  of  the  sth  wire,  for  s » 1,2 L 


(3.1b) 


i 
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d * 

where  E is  the  unknown  aperture  electric  field,  J.  is  the  surface 

inc 

current  density  on  the  2.th  wire,  is  the  /-component  of  the 

incident  magnetic  field  and  n * /y/e  is  the  characteristic  impedance 
of  the  space.  The  functions  S^,  S^,  and  are  defined  by 

00  £ 

S^(y,y’)  » I sink  V sink  y'  (3.2a) 

m,n»o  K*  -k*  ^ ^ 

mn 

. « k 

S,(x,y,y')  “ -5-  y — sink  x sink  y sink  y'  (3.2b) 

2 m,n.l  -k^  ^ y 

mn 

and 


S3(x,x',y,y')  = I 


m,n»l  K -k 
mn 


sink  X sink  x'  sink  y sink  y' 
XX  T T 


(3.2c) 


where 


k 


X 


k 


y 


nir 

b 


Note  that  any  of  the  double  sums  in  (3.2)  can  be  reduced  to  a single 
sum  by  use  of  (2.8)  and  will  thus  converge  exponentially  as  long  as 
f f'  [that  is,  when  (x,y)  ^ (x',y')]. 

Now  assume  that  no  two  wires  touch  one  another  and  that  no 
wire  is  in  the  aperture  region.  This  means  then  that  with  the  e.xception 
of  the  first  sum  all  of  the  infinite  sums  in  (3.1a)  are  uniformly  con- 
vergent. Indeed  this  first  sum  is  uniformly  convergent  everywhere 

rl 

except  when  y ■ y'  in  which  case  it  diverges  like  i—  . Similarly,  all 
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Slims  in  (3.1b)  are  uniformly  convergent,  with  the  exception  of  the 
term  where  Jl»s,  in  which  case  one  point  on  the  surface  of  integration 
will  coincide  with  the  point  (x,y). 


It  is  now  useful  to  make  the  approximation  that  the  wire  cur- 
rent resides  at  the  center  of  the  wire  and  that  field  boundary  condi- 
tions are  still  enforced  at  the  wire  surface.  This  approximation  is 
commonly  known  as  the  reduced  kernel  approximation.  For  this  problem 
it  can  be  stated  by 

S2(x,x',y,y')dc'  ~ 2TTr  S^Cx.x^.y.y^)  (3.3) 

■'wire 

where  (x^,y^)  is  the  center  point  of  the  wire  and  r is  the  wire  radius. 
This  can  be  justified  by  noting  that  S^,  the  Green's  function  for  the 
.magnetic  vector  potential  in  the  cavity,  must  contain  the  free  space 
two-dimensional  Green's  function,  which  goes  as  £n|r-r'|.  Note  that 
J!.n|r-r'|  can  be  integrated  analytically  over  the  wire  surface,  where 
r is  also  on  the  wire  surface,  to  give  the  desired  result  that 

rV 

2.n(2rsin  4^rdd  ' ■ 27rriln  r . 

Since  the  wire  is  thin  and  the  remaining  portion  of  S_  is  a smooth, 
homogeneous  solution  to  the  wave  equation,  it  too  can  be  validly  approx- 
imated by  this  technique.  It  should  be  noted  that  given  the  wire  radii 
small  but  finite,  the  application  of  the  reduced  kernel  approximation 
makes  all  of  the  infinite  sums  in  (3.1b)  uniformly  convergent. 
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Reduction  to  Matrix  Equation  by  Method  of  Moments 


An  effective  technique  for  the  numerical  solution  of  integral 


I 

equations  in  electromagnetics  is  the  method  of  moments  (Harrington 
1968).  Consider  the  operator  equation 


Lu  = F 


(3.4)  ‘ 


where  L is  a linear  operator,  F is  a known  vector,  and  u is  an  unknown 

1 

vector  for  which  the  solution  is  desired.  Using  the  method  of  moments,  ^ 

one  approximates  u by  a finite  linear  combination  of  expansion  vectors  | 

u . Thus  let  i 

n 


u 


N 


I 

q»i 


a u 


q q 


(3.5) 


Now  take  the  inner  product  of  (3.4)  with  N testing  vectors  w^.  If 
(3.S)  is  substituted  into  this  result,  and  it  is  noted  that  L is 
linear,  one  arrives  at  the  matrix  equation 


N 

[ <w  ,Lu  >a  = « P * 1,2,3,...N  . (3.6) 

q*i  p q q P 

• Note  that  (3.1a)  is  a linear  operator  equation  with  an  integro- 

differential  operator  and  unknown  function  E^(y)  and  the  unknowns 
Further,  one  has  the  boundary  condition  that  must  vanish  at  each 
edge  of  the  aperture  because  the  tangential  electric  field  (E^)  must 
vanish  at  these  two  points.  If  the  inner  product  for  this  space  is 
defined  to  be 
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fh+d 


<u,v> 


u(y)v(y)dy 


the  method  of  moments  may  be  applied  to  approximate  (3.1a)  by  a matrix 
equation. 

The  optimum  choice  of  expansion  functions  (u^)  and  testing 
functions  (w^)  for  this  integro-differential  operator  and  these  bound- 
ary conditions  has  been  the  subject  of  many  examinations  (Butler  and 
Wilton  1975;  Wilton  and  Butler  1976).  It  has  been  shown  that  one  effi- 
cient choice  (Wilton  and  Butler  1976)  is  that  of  pulse,  or  piecewise- 
constant,  expansion  functions  and  piecewise-sinusoidal  testing 
fiinctions,  denoted  p^  and  respectively  and  defined  in  the  coordinate 
system  of  this  problem  by 


Pq(y)  - N 


1 . |y-yqi  < I 


0.  ly-yJ>j 


(3.7) 


and 


Aj(y) 


sink(A-ly-ypl)  , [y-y^l  < A 


|y-ypl  > A 


(3.8) 


where  A » d/(N+l)  and  y^  » h+qA.  Thus  the  unknown  aperture  field  is 
approximated  by 


E^(y]  > I E P (y) 
q-1  ^ ^ 


(3.9) 
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Figure  6 shows  a plot  of  such  an  approximation  as  well  as  the 
testing  functions  A^Cy).  It  should  be  noted  that  a half-pulse  of  zero 
amplitude  has  been  placed  at  each  end  of  the  aperture.  This  is  done 
because  of  the  boundary  condition  that  E®  vanish  at  each  end  of  the 
aperture . 

In  order  to  perform  the  inner  product  one  takes  advantage  of 

the  piecewise-sinusoidal  testing  functions  and  the  differential  portion 

of  the  operator,  ^ k*,  and  integrates  by  parts  twice.  The  integral 
dy^ 

portion  of  the  result  vanishes,  leaving  only  three  boundary  terms. 

If  these  expansion  and  test  functions  are  applied  to  (3.1a)  and 
the  same  expansion  functions  are  substituted  into  (3.1b),  the  following 
equations  are  obtained 


X 


2CHj"‘=(o,yp) 


(3.10a) 


for  p “ 1,2, ... ,N  and 


N 

q-i  J 


S2(Xs,ys.y')dy'  - I IjS,(Xg,Xj,y^,y,)  = 0 


i-1 


rz'-  V 


(3.10b) 


3k 


for  s * 1,2,...,L  where 


Fq(y) 


Hf*^(k|y-y' |)dy'  ♦ j 


Aq 


Si(y,y')dy'  . (3.11) 


Aq 


In  (3.10),  * 2irrj^J^  is  the  current  due  to  the  uniform  surface  cur- 

rent density  on  the  2.th  wire  and  Aq  is  the  interval  (yq" 
point  centerpoint  of  the  £th  wire  and  = (x^.y^) 

is  a point  on  the  surface  of  the  sth  wire.  However,  for  thin  wires, 
this  distinction  is  important  only  when  Jl«s.  The  constant 


y +A 

:P 


C ■ 


Ap(y)‘ly  * 2(l-coskA) 


results  from  assuming  that  both  the  incident  magnetic  field  at  the 

aperture  and  the  magnetic  field  in  the  aperture  excited  by  the  currents 

on  the  wires  can  be  approximated  over  the  range  of  pth  testing  function 

by  their  respective  values  at  y^. 

It  is  now  observed  that  (3.10)  represents  two  coupled  matrix 

equations  for  the  unknowns  (E  } and  {ij}.  These  equations  can  be 

q X 

solved  simultaneously  by  standard  matrix  techniques,  such  as  Gauss 
elimination,  to  find  the  solution  for  the  unknown  aperture  field  and 
the  currents  excited  on  the  interior  wires. 

Numerical  Considerations 

Although  the  problem  has  been  formulated  and  the  integral 
equations  have  been  approximated  by  matrix  equation  (3.10),  it  is  noted 
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that  one  must  be  able  to  calculate  the  infinite  sums  S2  and  S^,  possibly 
near  the  source  where  convergence  is  poor.  The  sums  and  $2  must  be 
integrated,  over  the  source  region  (3.11).  Consequently,  it  becomes 
necessary  to  develop  techniques  for  handling  these  situations. 

First  it  should  be  noted  that  in  (3.10b)  it  is  valid  to  inte- 
grate $2  term  by  term  since  S2  will  always  be  uniformly  convergent 
over  the  range  of  the  integration  (Titchmarsh  1939,  p.  36).  It  can  be 
shown  that  in  (3.11)  can  also  be  integrated  term  by  term,  even 
though  at  y»y',  diverges.  Because  is  uniformly  convergent  at 
every  other  point  in  the  range  of  integration,  and  because  the  series 
that  results  if  Sj^  is  integrated  term  by  term  is  absolutely  convergent, 
it  is  valid  to  integrate  term  by  term  (Titchmarsh  1939,  pp.  44-4S) . 
This  is  an  important  result  since  as  a general  rule  it  is  numerically 
more  efficient  to  integrate  the  series  term  by  term. 

It  remains  to  develop  a technique  for  efficiently  Slimming  these 
series,  even  when  |r-r' | is  small.  It  is  expected  that  the  convergence 
will  be  poor  from  a numerical  standpoint.  To  demonstrate  the  method 
which  is  used  to  make  this  improvement,  consider  the  special  case  of 
Sj  (3.2c).  By  applying  (2.8a)  to  (3.2c)  one  finds  that 


2 “ sinhY^^y^sinhY^^(b-y^) 

S-  » — y sink  X sink  x'  v k 

3 a X X Y.  sinhr.b 

m»l  D b 


(3.12) 


where  » k^  - k^ . 

b X 
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Note  that  just  as  easily  the  double  sum  could  have  been  reduced 
by  eliminating  the  sum  over  m.  Thus,  there  is  always  a choice  of  which 
sum  to  remove  analytically  by  (2-8).  Let  be  defined  as  the  asymp- 
totic series  of  such  that  its  terms  are  the  limit  of  the  terms  in 
(3.12)  for  large  m.  Thus 


1 r 1 X . 

— ) — sinnnr—  sin 

*•,  m a 

m»l 


imrx' 


(3.13) 


Series  has  exponential  convergence,  the  rapidity  of  which  increases 
as  |y-y' I /a  increases.  It  can  easily  be  seen  that  if  the  siom  over  m 
were  removed  analytically,  would  converge  exponentially  like 
|x-x' 1/b.  It  is  reasonable  to  assume  from  this  that  as  a general  rule 
!y-y' |/a  should  be  compared  to  jx-x' |/b  and  (3.2)  should  be  reduced  to 
a single  sum  in  that  way  which  maximizes  this  exponential  convergence. 
For  example,  if  |y-y' |/a  > |x-x'l/b,  then  (3.2c)  should  be  reduced  to 
(3.12).  Indeed,  the  validity  of  this  general  rule  can  be  substantiated 
by  the  results  of  numerical  testing. 

Now  return  to  (3.12).  Suppose,  however,  that  although  |y-y' [/a 
is  larger  than  |x-x'|/b,  it  too  is  small,  and  (3.12)  remains  poorly 
convergent.  Convergence  of  such  a series  can  often  be  greatly  improved 
if  a closed  form  expression  can  be  found  for  the  corresponding  asymp- 
totic series  (Lewin  1975).  This  means  that  ifS»Zs,s  -*-s'as 
n -*•  »,  and  S'  * ^ ■ h where  h is  a closed  form  expression,  then 

D » H (s^  - s^')  will  converge  more  rapidly  than  S.  Consequently,  S 
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can  be  evaluated  by  S ■ h + D,  Thus  if  (3.13)  can  be  evaluated  in 
closed  form,  the  convergence  of  (3.12)  could  be  improved. 

For  this  particular  sum,  Sj,  one  rewrites  (3.13)  as 

asv  1 r e*®“ 

^3  * 2?  ^ “5“  (cosmg^  - cosme2)  (3.14) 

m“l 


where 


I 


ih-y' 


. h- 6,  ■ « 


But  it  is  known  that  (Jolley  1961,  pp.  110-111) 


-nx 


I ^ ^ cosnA  * j - y ^n(coshx  - cosA)  - y £n2  . (3. IS) 


n=l 


With  the  use  of  (3.14)  and  (3.15),  (3.12)  becomes 


2 " 

S_  a — / sink  X sink  x' 
3 a X X 

m"! 


‘sinhY^jy<sinhYj^(b-y^) 


Yb  sinhYbb 


1_  ^ 
2 k 


-ma 


j /cosha  - cos82\ 

4 IT  yco^a  - cos3^  j 


(3.16) 


It  should  be  noted  that  ntimerically  the  sum  in  (3.16)  is 
rapidly  convergent  regardless  of  how  small  |r-r' | is,  so  long  as  it 
is  not  zero.  Indeed  at  |r-r' | * 0,  a and  3^  go  to  zero  and  from  (3 
it  is  seen  that  possesses  the  expected  logarithmic  singularity. 


16) 
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Although  this  procedure  is  demonstrated  here  only  for  S^,  it 
can  also  be  successfully  applied  to  any  sum  or  term-by-term  integration 
of  a sum  which  is  needed  in  this  problem.  One  will,  however,  need  the 
following  asymptotic  series  in  closed  form: 


(3.17a) 


00 

^ sinnA 

n»l 


sinX 

2(cosKx  - cosX) 


(3. 17b) 


(3. 17c) 


(3.17d) 


where  is  the  ith  Bernoulli  number  as  defined  by  Abramowitz  and 
Stegun  (1965,  p.  804).  In  the  literature,  (3.17a)  can  be  found  in 

I 

Jolley  (1961,  p.  110-111),  (3.17b)  and  (3.17c)  are  found  in  IVheelon 

i 

(1968,  p.  38),  and  finally,  (3.17d)  is  found  in  Lewin  (1958,  p.  246) 
or  Lindelof  (1947,  p.  140).  In  (3.17d),  the  left  side  converges  well 
for  large  x and  the  right  side  converges  well  for  small  x. 


)le  Calculations 


Now  that  the  matrix  equation  has  been  formulated  and  methods 


of  computing  its  elements  have  been  devised,  solutions  for  various 

cavity  and  wire  configurations  can  be  obtain  via  numerical  solution 
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on  a digital  computer.  In  this  section,  a few  representative  solutions 
are  presented. 


It  should  be  noted  that  a similar  problem  in  which  the  infinite 
screen  is  omitted  can  be  solved  using  a formulation  and  computer  code 
previously  developed  (Seidel  1974)  to  calculate  the  currents  on  an 
array  of  cylindrical  scatters  in  free  space.  From  these  currents,  one 
can  easily  calculate  the  electric  field  in  the  aperture.  Therefore, 
in  addition  to  solutions  of  the  problem  at  hand,  solutions  to  the 
similar  problem  are  presented.  These  problems  are  referred  to  as  the 
flanged  and  unflanged  solutions,  respectively,  flange  meaning  the 
infinite  planar  screen  of  the  initial  problem. 

First  consider  the  case  of  a cavity  with  a depth  (a)  of  .6X 
and  a width  (b)  of  .8X.  This  cavity  has  an  aperture  width  of  .6A  which 
is  centered  in  the  cavity  wall  (d  » .6X,  h » .IX).  The  cavity  is 
excited  by  a plane  wave  with  unit  magnitude  electric  field  which 
impinges  on  the  cavity  from  the  negative  x direction.  Figure  7 shows 
a plot  of  aperture  field  for  this  cavity  for  the  case  of  no  internal 
wires.  It  is  seen  that  the  difference  between  the  flanged  solution 
and  the  unflanged  solution  is  relatively  small.  Since  the  interior 
fields  are  uniquely  determined  by  the  aperture  fields,  this  indicates 
that  the  presence  of  the  flanges  has  little  effect  upon  these  interior 
fields.  Note  also  that  the  aperture  fields  go  to  :ero  at  the  edges  of 
the  aperture  and  that  they  have  a maximum  magnitude  of  slightly  less 
than  unity.  The  fields  are  symmetric  about  the  center  of  the  aperture 
because  of  the  symmetry  of  the  cavity  itself  and  the  symmetry  of  the 
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incident  field. 


Consider  now  the  same  cavity  and  excitation,  but  with  a wire 
of  radius  .OOlX  located  at  the  geometric  center  of  the  cavity,  (.3X, 
.4X).  Figure  8 shows  the  aperture  field  for  this  case.  Again  the 
difference  between  unflanged  and  flanged  solutions  is  relatively  small. 
Note  that  this  is  also  true  for  the  currents  excited  on  the  wire,  which 
are  (-1.00  - j 1.48)ma  for  the  unflanged  case  and  (-.94  - j 1.38)ma  for 
the  flanged  case.  The  one  significant  difference  between  this  solution 
and  that  for  the  case  of  no  wires  is  that  the  maximum  magnitude  of  the 
aperture  field  is  approximately  three  times  larger  with  the  wire  than 
without. 

Again  consider  the  same  cavity  and  excitation,  but  now  with 
two  wires  of  equal  radii  (r  » .OOlX)  which  are  symmetrically  located 
about  the  center  of  the  cavity  at  (.3X,  .2SX)  and  (.3X,  .S5X). 

Figure  9 gives  the  aperture  field  for  this  case.  The  currents  on  the 
two  wires  are  equal  due  to  the  symmetry  and  are  (-.42  - j 1.19)ma  and 
(-.42  - j 1.2S)ma  for  the  unflanged  and  flanged  cases,  respectively. 

The  same  observations  that  were  made  for  the  case  of  one  wire  are  also 
applicable  here,  except  that  the  peak  magnitude  of  the  aperture  field 
is  for  this  case  slightly  larger. 

Finally  consider  the  same  cavity  and  excitation  for  one 
centrally  located  wire  (r  = .001),  but  with  a centered  aperture  of 
width  .4X,  The  aperture  fields  for  this  case  are  shown  in  Figure  10. 
The  wire  currents  are  (.82  - j 1.71)ma  and  (.98  - j 1.48)ma  for  the 
unflanged  and  flanged  cases  respectively.  By  comparison  to  Figure  8, 
it  is  seen  that  as  a result  of  shortening  the  aperture,  the  magnitudes 
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Aperture  Field  for  Cavity  with  Two  Interior  Wires. 
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of  the  aperture  field  and  current  have  slightly  Increased  and  the 
phase  has  changed  considerably.  Again  note  that  the  difference 
between  unflanged  and  flanged  solutions  is  relatively  small. 


Final  Observation^ 

At  this  point,  a summary  of  the  portions  of  this  two-dimensional 
problem  which  will  provide  insight  toward  the  three-dimensional  cavity 
problem  is  in  order.  Such  items  can  be  divided  into  two  categories: 
those  which  provide  insight  toward  efficient  computational  methods  and 
those  which  lead  to  a better  physical  understanding  of  the  problem. 

First  consider  those  items  in  the  computational  category.  If 
the  improvement  of  convergence  of  the  two-dimensional  sums  is  necessary 
in  this  problem,  then  surely  it  is  necessary  for  the  triple  sums  of  the 
three-dimensional  problem.  Such  techniques  as  removal  of  the  asymp- 
totic series  should  be  considered,  although  the  amount  of  work  done  in 
the  literature  on  double  sums  is  very  small  in  comparison  to  that  done 
on  single  series.  Another  numerical  technique  which  should  be  consid-  j 

ered  is  the  use  of  the  reduced  kernel  approximation  for  the  thin  wire 
scatterer.  Remember  that  this  approximation  removed  the  necessity  of 
performing  an  integration  of  a diverging  series. 

Secondly,  physical  insight  into  the  effects  of  modeling 
approximations  can  be  gained  by  noting  the  one  overwhelming  result  of 
the  data  presented.  The  interior  cavity  fields  and  the  wire  currents 
are  relatively  insensitive  to  the  presence  of  the  infinite  screen  in 
this  two-dimensional  analog  which  is  particularly  true  near  resonance 
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regions.  It  is  quite  reasonable  to  assume  that  this  will  also  be 
true  for  the  three-dimensional  problem,  and  thus  provides  valuable 
information  relating  to  modeling  the  three-dimensional  cavity 
problem  with  or  without  the  infinite  screen.  This  finding  is 
important  because  it  is  much  easier  to  consider  the  problem  with 
an  infinite  screen  than  without  a screen. 
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CHAPTER  4 


NUMERICAL  METHODS  FOR  CAVITY  PROBLEM 

Approximation  by  Matrix  Equation 
In  order  to  obtain  a solution  to  the  integral  equation  (2.13) 
for  the  three-dimensional  cavity  (with  the  coupled  constraints  in 
(2,21)),  it  is  possible  to  approximate  the  integral  equation  numeri- 
cally by  a matrix  equation.  This  is  generally  accomplished  by  the 
application  of  the  method  of  moments,  which  was  outlined  in  Chapter  3. 

Consider  the  integro-differential  equation  (2,13)  for  the 
unknown  current  J^(z).  It  is  seen  that  the  differential  portion  of 

the  operator  is  the  harmonic  operator  (^ — + k*),  which  was  encountered 

dz^ 

in  (3,1a).  Since  by  thin  wire  approximations  one  also  knows  that  the 
current  must  vanish  at  z and  Zn  (the  z-coordinates  of  the  wire  end- 

U *" 

points,  with  z > Zj),  the  use  of  piecewise  sinusoidal  testing  fume- 
tions  and  pulse  expansion  functions  is  indicated  (Wilton  and  Butler 
1976),  as  was  the  case  in  the  two-dimensional  problem.  For  this 
problem,  these  functions  are  defined  by 

ll  , |z-z  1 < 4 

I H - 

Pq(z)  * (4.1) 

jo  , |z-z  I > T 
\ q 

and 
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(4.2) 


sink(A-(z-Zp|)  , |z-Zp|  < A 

0 , |z-Zp|  > A 

where  A » (z^-Zj^)/ (N+1)  and  z^  » z^  + qA. 

Figure  11  illustrates  these  expansion  and  test  functions.  It 
is  important  to  note  that  if  the  wire  is  attached  to  the  cavity  at  one 
or  both  ends,  this  choice  must  be  slightly  modified  because  the  wire 
current  does  not  necessarily  vanish  at  the  attached  end.  To  rectify 
this  problem,  the  zero  half-pulse  at  an  attached  end  is  replaced  by  a 
half-pulse  of  unknown  amplitude.  Consequently,  a new  testing  function, 
which  is  a half-piecewise  sinusoidal  function,  must  be  introduced. 

These  are  shown  in  Figure  11  by  dashed  lines.  In  the  subsequent 
development  of  the  matrix  equation,  it  is  assumed  that  the  wire  is 
unattached.  However,  the  extension  to  attached  wires  is  straightfor- 
ward. 

Since  is  assumed  to  be  uniform  about  the  circumference  of 
the  wire  it  is  helpful  to  define  the  current  to  be  the  integral  of 
the  current  density  about  circumference  of  the  wire.  By  this  assump- 
tion 

I^(z)  - 2TTr  J^(z)  . (4.3) 

Now  approximate  by 

N 

1^(2)  - I Ipj:)  ■ (4.4) 

q«l  ^ ^ 


A;(z) 
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If  this  is  substituted  into  (2.13)  and  (2.21)  and  the  inner  product 
defined  by 


<u,v>  * 


u(z)v(z)dz 


is  taken  of  (2.13)  with  the  functions  A^,  the  following  coupled  matrix 
equations  are  obtained: 


I 

q 


0 . 


(4,5a) 


p * 1 , 2 , . . . , N , and 


(4.Sb) 


p * 1,2,3,  These  functions  provide  the  solution  for  the  unknown 
current  amplitudes  (I  ) and  the  unknown  aperture  fields  (t.,t-,t-)  » 

A A 0 

(E._  ,a-  ,H_  ) driven  by  the  exterior  short  circuit  fields  (E,  ,E-,E-)  » 
n ‘l  *2  1-3 

matrices  in  (4.5)  are  defined  by 


Q^q--jn[A^(Zp^P  - 2coskAA^(Zp)  > A^CJp.^)]  , (4.6) 

zn  ^ 

^2  ’ 
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(4.7c) 


0 - * - jCkna  (r„,r  ) 

^3  ”22  ®z2  P ^ 


and 


Q2q  - ’ 


1 ^ «eG^  - jkna^  g^  - jkna„  g 


nn 


®11  ®nl 


n — ® ^ 1 - Ct  g. 


ir\K 


®22  ®n2 


- ct  g. 

“22  ^12 


- a g,  1 - a 3. 

2n  ”11  ^21  ”22  ^22 


where  the  deleted  Green's  functions  in  (4.9)  are  evaluated 


» (r,,r  ).  Define 

SL  3 


A (z  ) 

q p 


K(Zp,2’)dz' 


Aq 


where  Aq  is  the  interval 


and 


l_ 

2:t 


J 

-IT 


G (?  .?')d(t)' 
zz  ^ 


(4.8a) 

(4.8b) 

(4.8c) 

> 

(4.9) 

at  (r,r') 

(4.10) 

(4.11) 


53 


where  (j»'  is  the  angular  coordinate  of  a cylindrical  coordinate  system 


about  the  wire  axis  (Figure  1).  Also,  C = 2(l-coskA)  and  Tp  and  r' 
are  on  the  wire  surface.  Note  that  to  obtain  (4.8),  the  integral  over 
the  qth  tubular  surface  segment  on  the  wire  has  been  approximated  by 
the  product  of  the  surface  area  of  the  segment  and  the  integrand 
evaluated  at  a point  on  the  center  of  the  segment. 

In  order  to  solve  the  coupled  matrix  equations  (4.5),  one  need 
only  obtain  a numerical  solution  to  the  partitioned  matrix  equation 


Note  that  if  the  effects  upon  the  aperture  field  of  the  fields  scat- 
tered by  the  wires  are  ignored,  = 5.  On  the  other  hand,  if  the 

effects  of  the  cavity  wall  reflections  upon  the  aperture  field  are 
sd  5 

Ignored  Q =1,  the  identity  matrix. 

However,  as  in  the  two-dimensional  problem,  difficulty  arises 

in  attempting  to  compute  the  elements  of  the  matrix,  each  of  which 

contains  one  of  the  triply-infinite  sums  defined  in  Tables  1 and  2. 

The  computations  which  exhibit  this  difficulty  can  be  categorized  into 

three  basic  types,  the  first  of  which  is  computation  of  any  one  of  the 

sums  outside  the  source  region.  The  second  category  is  that  of  the 

integral  of  G.  over  the  surface  of  the  wire,  for  which  at  one  point 
^22 

of  the  integration  G.  will  be  divergent.  Finally,  the  deleted 

*22 

Green's  functions  must  be  computed  at  the  aperture.  The  remainder  of 
this  chapter  attempts  to  deal  with  precisely  these  difficulties. 
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An  Efficient  Method  of  Computing  the  Sums 
Outside  the  Source  Region 

As  was  noted  in  Chapter  2,  each  of  the  Green's  functions 
(Tables  1 and  2)  can  be  reduced  from  a triple  sum  to  a double  sinn 
which  is  exponentially  convergent  for  Ir-r' | ^ 0 using  (2.8).  Indeed, 
it  is  easily  shown  that  the  asymptotic  series  associated  with  any  one 
of  these  exponentally  convergent  series  is  of  the  form 

-k  |z-2'  I 

= I f(m.n)  ^ 5 (4.13) 

m,n 

where  » (— ) + (^)  , a = 0,1,2  and  f(m,n)  is  a non-exponential 

function  of  rfi  and  n. 

As  was  found  in  Chapter  3,  it  is  a good  general  rule  to  reduce 
the  triple  sum  in  such  a way  as  to  produce  the  double  sum  with  the  most 
rapid  exponential  convergence.  For  example,  if 

(!_  * i_)(2.z.)2 

a^  b^ 

is  greater  than  both 

(—  + — )(x-x')"  and  (i-  ♦ l-)(y.y')^ 
b^  c^  a^  c^ 

then  the  sum  over  i should  be  reduced. 


Following  the  lead  from  the  two-dimensional  problem  in 

Chapter  3,  one  would  now  attempt  to  find  the  closed  forms  of  (4.13) 
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for  the  various  f(m,n)  and  a.  It  turns  out,  however,  that  the  list  • 

of  double  sums  for  which  a closed  form  is  known  is  painfully  short, 
and  does  not  seem  to  include  (4.13)  for  any  a and  f(m,n)  of  interest 
in  this  problem. 

Another  method  which  has  been  suggested  is  to  remove  the 
known  singularity  from  the  series  by  expanding  the  singularity  in  the 
same  expansion  functions  as  the  series  itself.  Both  Tai  and  Rozenfeld 
(1976)  and  Rahmat-Samii  (1975)  have  removed  a delta  function  singular- 
ity from  their  series  for  5^  in  this  way.  However,  this  serves  little 
purpose  because  the  kpown  singularity  of  is  not  the  delta  func^on. 

In  fact,  as  shown  by  Howard  (1974),  the  singularity  of  5^  is  actually 
the  longitudinal  portion  of  l6(f-f'). 

The  third  alternative,  and  the  one  that  is  used  for  this  work, 
is  to  simply  sum  the  double  series  in  an  efficient  manner.  Note  that 
because  of  the  exponential  convergence,  as  one  attempts  to  make  this 
computation  nearer  and  nearer  the  source,  the  series  will  become  more 
and  more  poorly  convergent.  Thus  one  should  expect  to  reach  a point 
such  that  for  |r-f ' | less  than  some  minimum  value,  numerical  computa- 
tion of  the  sum  in  this  fashion  becomes  unfeasible. 

However,  some  things  can  be  done  which  make  this  method  more 
efficient.  Because  of  the  exponential  convergence  in  the  asymptotic 
series  (4.13),  one  would  expect  an  efficient  ordering  of  terms  to  be 
in  order  of  increasing  k^.  This  takes  advantage  of  the  exponential 
convergence  as  well  as  the  k^  in  the  denominator.  At  this  point,  it 
is  useful  to  partition  the  m-n  plane  with  successive  curves 
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(Bromwich,  1926,  p.83).  If  the  sum  of  all  terms  lying  between  two 


successive  curves  is  called  s^,  then  the  double  series  can  be  converted 
into  the  single  series  of  the  form 


oa 

I s 

q-1  “ 


By  the  proper  selection  of  these  curves,  the  most  efficient  ordering 
of  terms  can  be  determined. 

For  this  problem  such  a choice  would  be  that  of  ellipses  with 
semi-axes  in  m of  a a and  semi-axes  in  n of  a b,  where  a is  a monoton- 

<i  q q 

ically  increasing  sequence  of  constants.  Such  a partitioning  is  shown 
in  Figure  12.  Note  that  for  such  a choice,  each  successive  partition 
contains  terms  for  which  is  larger  than  in  the  preceding  partition. 
Also  note  that  since  the  sum  of  terms  in  the  qth  partition  is  the  qth 
term  of  a single  infinite  series,  methods  used  for  determining  the 
convergence  of  single  series  can  be  applied. 

To  test  this  method,  the  sums  were  numerically  computed  via 
digital  computer.  Figure  13  shows  the  notation  used  for  the  sum. 

Note  that  it  is  assumed  without  loss  of  generality  that  a ^ b. 
Numerically  the  series  was  truncated  to  include  only  those  terras  with- 
in C^.  The  maximum  value  of  m was  M.  Note  that  the  total  number  of 
terms  is  approximately  the  area  within  C,^  or  N.^.  ~ ^ . Let  » 

P 00 

T s and  R * T s . Then  R,,  is  the  error  resulting  from  trunca- 
1 q P 1 q •' 

tion.  R^  can  be  crudely  bonded  from  above  with  an  integral  bound. 
However,  for  actual  computation,  convergence  was  defined  to  have  been 
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reached  when  ratio  (s  /S  ) was  less  than  some  small  constant  C for 

P P s 

Z consecutive  values  of  o. 


s 

Note  that  any  one  of  the  Green's  functions  can  be  reduced  to 
three  entirely  different  double  sums.  For  the  purpose  of  testing,  the 
Green's  functions  were  computed  by  all  three  double  sums,  and  then 
these  results  were  compared.  Indeed,  the  effectiveness  of  this  test 
is  attested  to  by  the  several  programming  errors  which  it  detected. 
However,  once  these  errors  were  corrected,  extensive  testing  demon- 
strated that  the  three  values  always  agreed  to  appro.ximately  the 
accuracy  specified  by  the  constant  C^.  This  test  also  showed  that  the 
general  rule  of  reducing  the  sum  so  as  to  give  the  most  rapid  exponen- 
tial convergence  did,  in  fact,  produce  the  particular  one  of  the  three 
possible  double  sums  which  required  the  fewest  number  of  terms. 

Finally,  extensive  testing  demonstrated  that  for  cavity  sites 
in  the  vicinity  of  the  first  resonance,  the  Green's  functions  could  be 
easily  obtained  for  values  of  jr-r' | greater  than  A/20  (for  larger 
cavities,  this  minimum  distance  increases;  for  smaller  cavities,  it 
decreases).  Computations  at  even  smaller  values  of  |r-f'|  are  not 
impossible,  but  rather  more  and  more  time  consuming. 


Numerical  Evaluation  of  the  Integral 
of  the  Singular  Sum 

Consider  now  that  evaluation  of  defined  by  (4.10)  and 

(4.11).  As  earlier  noted,  when  p « q,  the  integrand  G^  of  the 

zz 

integral  over  the  tubular  wire  surface  segment  diverges  at  = r'. 

Even  for  p ^ q,  if  p is  near  q then  the  integrand  will  converge  poorly. 
These  two  difficulties  must  be  overcome. 
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In  order  to  sidestep  the  first  problem,  it  would  be  helpful 


to  apply  the  reduced  kernel  approximation  (assume  all  current  resides 
at  center  of  wire).  With  this  approximation,  when  p = q the  integrand 
would  never  diverge  and,  in  fact,  would  be  uniformly  convergent  every- 
where on  the  surface  of  integration  because  |r  -r' | > r.  However,  one 
very  important  consideration  is  the  validity  of  using  the  reduced 
kernel  approximation  for  the  kernel  (4.11).  This  kernel  is  known  to 
include  the  singularity  of  the  free  space  kernel  plus  a remaining 
smooth  homogeneous  solution.  If  r « X,  it  is  justifiable  to  assume 
that  the  smooth  part  of  K is  essentially  the  same  at  the  center  of 
the  wire  and  at  points  on  *the  wire  surface.  Thus  one  needs  a like 
comparison  for  the  singular  portion. 

Define  K and  K to  be  the  free  space  exact  and  reduced 
or 

kernels,  respectively,  given  by 


K^(5) 


[5^  ♦ 4r^ 


sin 


*,-1/2 

2-* 


d* 


(4.14a) 


and 


Ky(5)  » [C"  > . 


(4.14b) 


Since  ultimately  integrals  of  the  kernels  over  5 are  needed,  such 
integrals  will  be  compared.  Let 


*,(2) 


K^(^)d5  . I^ 


6o 


(4.15a) 


and 


'i'j.Cz) 


'Z 

= ln[f  > 
• o 


C4,lSb) 


where  has  been  defined  by  Butler  (1975].  Figure  14  shows  a compari 
son  of  these  two  functions.  It  is  readily  apparent  that  such  an 
approximation  is  valid  for  the  singular  part  of  the  kernel  and  thus 
valid  for  the  kernel  (4.11). 

Using  the  equation  for  G.  from  Table  1 (reducing  it  to  a 

^zz 

double  sum  by  (2.8b))  and  applying  the  reduced  kernel  approximation  to 
(4.11),  one  obtains 


® coshY^z.coshy  (c-z  ) 


m,nal 


Y^sinhy^c 


F(x^,y^)  (4.16) 


where  , 

c X y 


F(x  ,y  ) = sink  X sink  x sink  y sink  y 
^ c’-'c'^  X c X p rc  y'p 


(4.17) 


(x^,y^)  is  the  location  of  the  center  of  the  wire  and  r^  = (.x^  + rcos<(), 

y + rsin<t,  z ) is  a point  on  the  wire  surface.  By  using  a hyperbolic 
c p 

trigonometric  identity,  one  can  express  (4.16)  as  the  sura  of  two  terras 
by 


K(Zp,z')  = S(|Zp-z' I)  + S(Zp>z') 


(4.18) 


where 


6l 


(4.19) 


S(8) 


2 


I 


cosh  Y (c-0) 

- sinh  Y c 
c c 


To  evaluate  A^C^p)  defined  by  (4.10),  one  needs  to  calculate  the 
integral  with  respect  to  z'  of  the  kernel  (4.18)  and  thus  of  (4.19). 
If  this  integration  is  performed  on  (4.19)  term  by  term,  and  sub- 
stituted into  (4.18)  and  (4.10),  one  finds  that 


A (z  ) » Qdz  -z  j)  Q(z  >z  ) 

q'-p'^  xvip  p q 


(4.20) 


where 


Q(a)  * ( 


P(B) 


p(e) 


A 

“ * 2 A 

« A ' “ -I 

8 ■ a - y 


I- 

6-0 


> P(6) 


2 - “ 
6-0 


' “ 1 I 


(4.21) 


and  where  P(6)  is  the  indefinite  integral  of  S(6)  given  by 


P(B) 


2 

iF 


00 


I 

m,n-l 


sinh  Y (c-6) 

^ ^ • 
Y*  sinh  Y c 
c c 


(4.22) 


Thus,  if  P(6)  can  be  evaluated  for  6^0,  then  Aq(Zp)  can  be  evaluated 
on  the  wire  using  the  reduced  kernel. 

Consider  the  case  where  6-0.  .Note  that  the  hyperbolic  sine 
functions  in  (4.22)  cancel,  leaving 
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P(0) 


(4.23) 


- <»  sink  X sink  x sink  y sink  y 

2 ^ X c x^p  fz  y^p 


— y 

ab  ^ 


m,n=l 


k^  + k^  - k^ 
X y 


I It  is  now  recognized  that 

I 

I 


P(0) 


- j S3Cx^,Xp,y^,yp) 


(4.24) 


where  is  defined  by  (3.2c).  Thus  P(0)  can  be  readily  evaluated 
numerically  using  (3.16)  for  r > o.  Note  that  in  (4.22),  because  the 
hyperbolic  sine  is  an  odd  function, 

PC2c  - 8)  - -P(8)  . (4.25) 

Thus  by  (4.24)  and  (4.25)  it  is  also  true  that 

P(2c)  - ^ S3(x^,Xp,y^,yp)  . (4.26) 

For  S>*0  or  2c,  there  appears  to  be  no  alternative  other  than  to 
perform  the  two-dimensional  sum  by  the  methods  of  the  previous  section. 
Note  that  asymptotically  P is  of  the  following  form: 

as  -k  8 

pasy  ^ I p e_c_  ^ 0 < 3 < ^ . 

m,n“l  k* 
c 

Thus  for  8 near  zero  (or  near  2c,  as  indicated  by  (4.25)),  poor  conver- 
gence is  expected. 

It  is  known  that  the  reduced  kernel  (4.16)  must  contain  the 

singular  portion  free-space  reduced  kernel  (4.14b)  plus  a smooth 

6U 


homogeneous  function.  By  (4,18),  S(g)  must  also  contain  that 
singularity.  If  a function  ij;(z)  is  defined  similar  to  (4.15b)  by 


'l'(z) 


•z 

S(8)de  =■  P(z)  - P(0) 
• o 


(4.27a) 


and  4*2(2)  is  defined  by 

4'^(z)  » - tli^(z)  (4.27b) 

then  4*2  is  the  integral  of  a smooth  function  and  thus  itself  is  smooth. 

Since  P(0)  is  readily  computed,  'J;(z)  can  be  computed  numeri- 
cally for  z greater  than  some  minimum  value  z . Note  that  the  function 

0 

ip  (z)  can  be  calculated  by  using  (4.1Sb)  for  any  z.  Thus  4*_(z)  can  be 

numerically  evaluated  for  z > z^.  Note  that  if  4*2  is  smooth  and  z^  is 

sufficiently  small,  4*2(2)  can  be  interpolated  for  o < z < z^.  Then  if 

’l*j.(z)  is  added  to  these  interpolated  values  of  ’^*2(2)*  4'(z)  can  be 

found  for  o < z < z . 

0 

To  demonstrate  the  practicality  of  this  technique,  consider  the 

curves  of  Figure  IS.  They  show  the  functions  ’Kz) , 4*^(2)  and  '4*2(2)  for 

o £ z £ .5X.  The  cavity  dimensions  are  a = .7X,  b =»  .8X  and  c » .3X 

with  a wire  of  radius  4 ■ .OOlX  located  at  x » .35X  and  y = .4X. 

c c 

Note  that  for  z < .2X,  is  almost  linear  and  could  be  interpolated 
quite  accurately.  The  utility  of  this  method  is  appreciated  if  it  is 
noted  that  for  this  example  the  calculation  of  P(z)  took  98  terms  at 
z » .2X,  242  terms  at  z » .IX  and  2139  terms  at  z * .02SX.  In  each 
case,  the  convergence  criteria  was  ■ 10  * and  « 3. 
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Before  leaving  this  section,  one  important  observation  should 

be  made.  By  examination  of  (4.20),  it  is  seen  that  A (z  ) is  a func- 

R P 

tion  of  two  terms,  one  of  which  depends  only  upon|p-q|  and  the  other  of 
which  depends  only  upon  (p+q) . This  means  that  A^(Zp)can  be  calculated 
for  all  values  of  q and  p by  3N  + 1 computations  rather  than  + 2N 
calculations.  Thus  for  N ^ 2,  the  computing  time  required  can  be 
significantly  decreased. 

Numerical  Evaluation  of  the  Deleted  Green's  Functions 
In  the  previous  sections  of  this  chapter  all  of  the  computation- 
al difficulties  in  filling  the  matrix  in  (4.12)  have  been  resolved  in  a 
workable  fashion  except  for  evaluation  of  the  deleted  Green's  functions 
in  (4.9).  These  functions,  defined  by  (2.17)  and  (2.18),  must  be 
evaluated  at  the  aperture  (r*r  ) . Note  that  although  these  deleted 
functions  are  bounded  solutions  to  the  homogeneous  wave  equation  at  the 
point  r , they  cannot  be  calculated  directly  from  (2.17)  because  both 

A 

2.(r,r  ) and  g_(r,r  ) are  divergent  at  . 

Ideally,  one  would  like  to  have  expansions  for  the  various 
free-space  dyadic  Green's  functions  so  that  the  singularities  could  be 
removed  from  the  cavity  dyads  term  by  term,  as  was  discussed  earlier 
in  this  chapter.  The  resulting  convergent  sums  would  be  precisely  the 
deleted  Green's  functions  of  (2.17)  and  (2.18).  However,  at  the  pre- 
sent time,  practical  techniques  of  employing  this  method  are  not 
available. 
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In  this  work,  a different  approach  was  taken  that  relies  upon 
some  of  the  earlier  observations  of  this  chapter.  Since  the  deleted  / 

Green's  functions  are  homogeneous  solutions  to  the  wave  equation,  it 
is  logical  to  expect  that  they  could  be  found  by  using  an  interpola- 
tion method  similar  to  that  employed  in  the  previous  section  to 

evaluate  A (z  }. 

q P 

Note  that  such  an  interpolation  scheme  would  require  the 
evaluation  of  several  different  components  of  the  various  cavity  dyads 
at  points  near  the  source.  If  it  is  assumed  that  the  accuracy  of  the 
interpolations  would  improve  as  the  points  move  nearer  the  source, 
increased  accuracy  would  require  increased  computation  time. 

Since  the  double  sums  in  the  dyads  for  the  vector  potentials 
converge  more  rapidly  than  those  for  the  fields,  it  is  advantageous  to 
compute  only  the  components  of  5^  and  gp  at  several  points  near  the 
source.  Then  by  (2.17),  the  associated  deleted  Green's  functions  for 
the  vector  potentials  could  be  subsequently  computed.  To  find  the 

/S  A 

values  of  5^  and  gp  in  the  aperture,  an  interpolation  scheme  was  em- 
ployed. To  find  the  values  of  the  other  deleted  Green's  dyads,  the 
differential  operators  in  (2.18)  were  approximated  by  finite  difference 
• techniques.  Note  that  the  points  at  which  the  dyads  are  computed 

should  be  chosen  carefully  so  that  they  provide  the  proper  information 
for  the  interpolation  and  finite  difference  techniques. 

A crucial  criterion  for  the  success  of  this  method  is  the 
smoothness  of  the  deleted  Green's  functions.  Figure  16  shows  plots  of 

A 

G,  and  Re(G^  ) as  a function  of  distance  from  the  aperture.  Note 
''zz  ''zz 
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that  the  deleted  function  is  indeed  smooth  and  could  be  approximated 
near  the  source  by  interpolation  techniques.  Figure  17  shows  similar 


plots  of  g and  ReCg-  ) and  the  same  observations  can  be  made. 

yy  yy 

This  method  has  been  tested  numerically  for  several  aperture 
location/ cavity  size  configurations.  For  each  configuration,  the 
computations  were  made  with  the  distances  between  the  aperture  and  the 
points  at  which  the  computations  were  made  becoming  successively 
smaller.  The  calculated  values  of  the  deleted  Green's  functions  for 
the  successively  smaller  distances  were  then  compared  in  order  to 
determine  if  the  process  was  convergent.  In  all  cases  tested,  it  was 
indeed  convergent.  Although  error  bounds  are  not  available,  based 
upon  the  numerical  testing,  accuracy  to  within  10  to  15  percent  is 
estimated  for  distances  from  aperture  to  computation  points  being 
approximately  .05X  to  .IX. 

It  should  be  noted  that  these  values  for  the  deleted  Green's 
functions  are  computed  with  less  accuracy  than  any  other  terms  in  the 
matrix  of  (4.12).  Nonetheless,  their  computation  requires  far  more 
computer  time  than  any  of  the  other  terms  in  the  matrix.  In  addition, 
because  these  functions  are  always  multiplied  by  the  aperture  polar- 
izabilities (which  are  small  for  small  apertures)  in  (4.9),  they 
manifest  themselves  primarily  as  perturbations  and  have  only  a slight 
effect  upon  the  solution,  as  will  be  seen  in  the  following  chapter. 
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Figure  17.  Example  Computation  of  g„  and  the  Real  Part  of  g, 


CHAPTER  S 


NUMERICAL  RESULTS  FOR  SAMPLE  CASES 

A computer  code  has  been  written  which  implements  the 
numerical  formulation  presented  in  the  previous  chapter.  In  this 
chapter,  selected  numerical  results  are  presented  which  demonstrate 

\ 

the  capabilities  of  this  numerical  code.  The  program  allows  for  an 
elliptic  aperture  in  any  one  of  the  three  walls  defined  by  x = 0, 

I 

y a 0 or  2 * 0,  and  assumes  that  the  wall  containing  the  aperture  is 

an  infinite  planar  screen.  For  all  cases  presented  here,  the  aper-  i 

ture  was  chosen  to  perforate  the  x = 0 wall  of  the  cavity.  Also  ' 

note  that  all  lengths  are  in  units  of  wavelength  A. 

I 

As  a first  example,  consider  a relatively  small  cavity 
whose  dimensions  are  defined  by  a * .2,  b = .25  and  c = .3  which  . 

encloses  a wire  of  length  .2  and  radius  .001.  The  wire  is  exactly  ^ 

centered  in  the  cavity  (z,  = .05,  z » .25,  x » .1,  y = .125). 

X(  U C C ‘ 

The  aperture  is  circular  with  a radius  of  .01  and  has  its  center  at  \ 

r » (0.,  .125,  .15).  The  exterior  excitation  is  a plane  wave 

a I 

1 

which  impinges  from  the  ^ direction  and  is  polarized  such  that  the  | 

I 

electric  field  is  in  the  ^ direction. 

Figure  18  shows  the  current  excited  upon  the  wire  for  this 

configuration.  Since  both  aperture  and  wire  are  symmetric  about  the 
c 

plane  2 “ j,  one  would  expect  to  observe  symmetric  properties  in  the 


72 


wire  currents.  This  is  indeed  the  case.  The  real  part  of  the 
current,  which  is  excited  by  the  magnetic  aperture  dipole  moment, 
possesses  even  symmetry,  where  the  imaginary  part,  which  is  excited 
by  P^,  has  odd  symmetry  about  z = .15.  This  relationship  between 
the  real  and  imaginary  parts  of  the  current  and  and  P^  respec- 
tively can  be  seen  by  examination  of  C4.Sa),  (4.6)  and  (4.7).  Also 
note  the  magnitude  of  the  current  which  peaks  at  approximately  .04  pa. 
One  would  expect  currents  on  the  order  of  1 ma  for  the  same  wire  in 
a free  space  environment.  Thus  the  shielding  of  the  cavity  reduces 
currents  by  approximately  five  orders  of  magnitude. 

Next  consider  the  case  of  a much  larger  cavity  containing  a 
longer  wire.  For  this  case,  a = .4,  b = .6  and  c = 1.3,  which  is 
larger  than  the  first  several  cavity  resonances.  Again  the  wire  is 
centered  in  the  cavity  =.15,  z^  = I.IS,  x^  = .2,  y^  = .3)  and  is 
one  wavelength  long  with  radius  r = .001.  The  elliptic  aperture  has 
semiaxes  of  .05  amd  ,01  in  the  y and  z directions,  respectively  and 
is  located  at  r = (0,  .2,  .4).  Again,  the  incident  plane  wave 
impinges  from  the  ^ direction  with  a directed  electric  field. 

Since  the  wire  is  of  resonant  length,  one  would  expect  to 
excite  resonant  currents.  Indeed,  as  shown  in  Figure  19,  this  is  the 
case.  Also,  note  that  the  current  magnitude  now  peaks  at  approxi- 
mately 10  pa.  This  increase  over  the  previous  case  can  be  attrib- 
uted to  t.hree  causes.  First,  the  aperture  is  larger  and  thus  more 
energy  is  coupled.  Second,  the  wire  is  of  resonant  length.  Finally, 
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5 


Figure  19.  Currents  Excited  on  lA  Wire  in  .4A  x .6A  x 1.3A  Cavity. 


the  smaller  cavity  is  considerably  below  the  lowest  cavity  resonance, 
whereas  the  larger  cavity  is  near  several  resonances. 

Note  that  if  one  semiaxis  of  the  aperture  is  much  larger  than 
the  other,  the  aperture  begins  to  look  like  a short  slot.  One  would 
expect  that  the  strongest  coupling  would  occur  when  the  slot  is  per- 
pendicular to  the  wire  and  the  incident  electric  field  is  perpencUcular 
to  the  slot.  To  test  this,  consider  a cavity  with  dimensions  .7  x .7 
X .8  (a,  b and  c respectively)  with  a one-half  wavelength  wire  of 
radius  r = .001  which  is  located  in  the  cavity  at  » .15,  = .65, 

X = .35,  y = .4.  The  aperture  is  located  in  the  x » 0 wall  [r  = 

(0.,  .4,  .4)],  and  has  semiaxes  of  .07  and  .01.  The  plane  wave  is 
incident  from  the  ^ direction  (normal  to  wall  of  aperture)  and  has  a 
£ directed  electric  field. 

Consider  two  cases:  that  where  the  slot  is  perpendicular  to 

the  incident  electric  field,  and  that  where  the  slot  is  parallel  to 

the  incident  electric  field;  that  is,  where  the  major  semiaxis  of  the 

aperture  is  in  the  y or  z direction,  respectively.  Figure  20  shows 

the  current  excited  upon  the  wire  for  these  two  cases.  It  is  readily 

•inc 

seen  that  the  current  excited  when  the  slot  is  perpendicular  to  E 
is  approximately  twenty  times  larger  than  that  excited  when  the  slot 

• 4 TIQ 

and  E are  parallel.  Thus  the  expected  effects  are  observed. 

Now  consider  the  same  case  except  change  the  polarization  of 
the  incident  plane  wave  such  that  the  electric  field  is  in  the  ^ 
direction.  Whereas  in  the  previous  case  the  incident  electric  field 
was  perpendicular  to  the  slot,  it  is  now  parallel  to  the  slot,  and 
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Fi^jure  20.  Currents  Excited  on  .5X  Wire  for  Slot  Perpendicular  and 
Parallel  to  : Directed  Incident  Electric  Field. 


consequently  one  would  expect  the  excited  currents  to  be  much  smaller. 
As  shown  in  Figure  21  this  is  the  case.  Again  the  currents  are  shown 
for  both  the  slot  perpendicular  and  parallel  to  the  incident  electric 
field.  Note  that  the  currents  are  again  related  by  a factor  of  ap- 
proximately twenty,  but  are  over  five  orders  of  magnitude  less  than  the 
currents  shown  in  Figure  20.  However,  for  the  same  wire  and  excitation 
in  free  space,  if  thin  wire  approximations  are  used,  there  will  be  no 
current  excited  upon  the  wire. 

An  examination  of  (4.5)  shows  the  small,  but  nonetheless 
nonzero  current  in  Figure  21  is  due  to  the  depolarizing  effects  of 
the  cavity  wall  reflections.  Indeed,  if  one  lets  = T in  (4.5), 
this  polarization  will  not  excite  the  wi^'e  in  the  cavity.  Since  for 
small  apertures,  these  effects  manifest  themselves  as  perturbations, 
one  does  expect  the  currents  excited  by  this  depolarizing  effect  to 
be  small. 

Another  effect  of  the  cavity  wall  reflections  is  seen  by 
noting  that  for  the  case  in  Figure  20,  since  the  excitation  is  nor- 
mally incident,  (r  ) is  zero,  and  thus  if  the  wall  reflections 
are  neglected,  would  be  zero.  However,  the  wall  reflections 
produce  a small  but  nonzero  and  thus  a small  imaginary  current. 
Although  not  shown  in  Figures  20  or  21,  there  was  actually  an  imagi- 
nary part  of  the  current,  which  in  all  cases  was  more  than  two  orders 
of  magnitude  less  than  the  corresponding  real  portion. 

At  this  point  it  is  helpful  to  ascertain  the  size  of  this 
perturbation  effect.  In  order  to  accomplish  this,  solutions  were 
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Figure  21.  Currents  Excited  on  .5X  Wire  for  Slot  Perpendicular  and 
Parallel  to  -y  Directed  Electric  Field. 


obtained  in  two  additional  manners.  In  the  first,  the  cavity  wall 
reflections  were  neglected  (5*^  = ?) . In  the  second,  the  wire  scatter 
ing  effects  at  the  aperture  were  also  ignored  (Q  = 0,  Q = I) . 

These  solutions  were  then  compared  for  several  cases,  including 
those  presented  in  this  chapter.  It  was  found  that  the  difference 
between  solutions  increased  as  the  aperture  size  increased.  This  is 
expected  as  the  perturbation  is  on  the  order  of  the  aperture  polariz- 
abilities, which  increase  as  the  size  of  the  aperture  increases. 
However,  in  no  case  did  the  difference  between  the  original  solution 
and  the  first  additional  solution  exceed  one  percent  (largest  aper- 
ture considered  was  circular  with  radius  .05)  and  the  difference 
was  much  less  for  smaller  apertures.  The  difference  between  the  two 
additional  solutions  was  even  smaller. 

It  is  apparent  that  for  these  small  apertures,  a solution 
closely  approximating  the  e.xact  one  can  be  obtained  by  simply 
neglecting  these  perturbation  effects.  Since  for  apertures  much 
larger  than  those  considered,  the  applicability  of  small  aperture 
theory  becomes  increasingly  questionable  anyway,  this  result  is 
quite  significant.  Since  the  computation  of  these  perturbations  are 
very  time-consuming,  their  omission  would  substantially  decrease 
computing  time. 

Finally,  consider  the  case  of  wires  connected  to  the  cavity 
at  one  or  both  ends.  Figures  22  and  23  show  the  wire  currents  for 
these  two  cases  respectively.  The  cavity  size  is  .7  x .8  x .8  and 
the  wire  axis  is  at  (.x^,  y^)  = (.IS,  .5)  with  a wire  radius  r = .001. 
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The  aperture  is  located  at  r = CO.,  .3,  .6)  and  its  semiaxes  in  the 

a 

y and  z directions  are  .07  and  .01  respectively.  The  plane  wave  is 
normally  incident  upon  the  aperture  from  the  -x  direction  and  the 
electric  field  is  z directed. 

In  Figure  22  the  wire  is  connected  at  z = .8  and  as  expected, 
the  axial  (z)  derivative  of  the  current  goes  to  zero  at  the  wall. 
Also,  as  expected,  the  current  at  the  free  end  of  the  wire  vanishes. 
Similarly  in  Figure  23,  where  the  wire  is  connected  at  both  ends  of 
the  cavity,  the  axial  derivative  of  the  current  vanishes  at  both 
z = 0 and  z = .8.  In  both  cases,  the  current  magnitudes  are  on 
the  order  of  10  ya.  For  a similar  cavity  and  excitation,  but  hav- 
ing a free  wire  (Figure  20,  perpendicular  slot)  the  current  magni- 
tude is  also  of  this  order. 

Thus  it  has  been  shown  that  physically  reasonable  numerical 
solutions  can  be  obtained  for  a variety  of  cavity/aperture/wire  con- 
figurations. It  should  be  noted  however  that  important  comparisons 
between  theory  and  experiment  must  await  the  publishing  of  experimen- 
tal results  applicable  to  this  problem.  Such  comparisons  would 
demonstrate  the  applicability  of  the  modeling  and  the  accuracy  of  the 
solution. 
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CHAPTER  6 


CONCLUSION 

In  this  paper,  the  task  was  undertaken  to  formulate  and  numeri- 
cally solve  the  problem  of  an  aperture  excited  wire  scatterer  in  a 
rectangular  cavity.  The  formulation  in  Chapter  2,  although  tedious, 
was  relatively  straightforward.  It  was  found  that  the  major  difficul- 
ties of  the  problem  lay  in  the  numerical  evaluation  of  the  infinite 
double  sum  Green's  functions  for  the  cavity  interior. 

By  considering  the  preliminary  two-dimensional  problem  in 
Chapter  3,  valuable  insight  was  obtained  toward  resolving  the  numerical 
difficulties  of  the  three-dimensional  problem.  It  is  important  to 
realize  that  the  two-dimensional  problem  is  also  a significant  problem 
in  and  of  itself.  Two-dimensional  problems  in  many  cases  provide 
adequate  models  for  more  complex  three-dimensional  structures.  In 
addition,  solutions  to  this  problem  have  the  advantage  of  being  appli- 
cable to  large  apert;ires  whereas  the  major  solution  herein  is  restricted 

to  small  apertures  only. 

The  real  significance  of  this  work  is  contained  in  Chapter  4. 
Here  it  was  demonstrated  that  the  dyadic  Green's  functions  for  the 
cavity  problem  can  be  calculated  (although  in  some  cases,  only  with 
considerable  effort).  This  finding  is  particularly  essential  in  or 
near  the  source  region.  This  is  because  the  treatment  of  the  singular 
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kernel  of  an  integral  equation  in  the  source  region  is  crucial  to  its 
numerical  solution.  Thus  the  ability  to  numerically  solve  the  integral 
equation  for  a scatterer  in  a cavity  is  demonstrated,  and  indeed, 
numerical  results  can  be  provided  as  found  in  Chapter  S. 

Suggested  Extensions  of  this  Work 
Perhaps  more  important  than  what  this  work  has  accomplished  is 
rather  what  extensions  and  applications  can  be  found  for  it.  In 
general,  the  most  immediate  extensions  of  this  work  would  be  to  elimi- 
nate some  of  the  restrictions  caused  by  the  initial  assumptions  of  the 
problem.  The  relevance  of  an  extension  can  thus  be  measured  by  asking 
how  restrictive  is  the  assumption  which  the  extension  eliminates. 

A very  important  extension  would  be  to  allow  a large  aperture 
in  the  cavity  wall.  It  should  be  noted  that  in  order  to  do  this  an 
aperture  field  integral  equation  must  be  formulated  and  solved  simul- 
taneously with  the  integral  equation  for  the  wire  (2.13).  In  general, 
the  aperture  fields  must  be  divided  into  a two-dimensional  array  of 
surface  patches,  requiring  large  amounts  of  computer  storage  and  time. 
However  if  the  aperture  is  small  in  one  dimension  (a  slot),  it  need 
only  be  divided  into  a one-dimensional  array  of  surface  patches. 

A second  extension  of  the  present  work  would  be  to  account  for 
more  complex  scattering  geometries.  This  might  include  uniformly  or 
lumped  loaded  wires,  more  than  one  wire,  or  wires  which  are  not  paral- 
lel to  one  of  the  coordinate  axes.  The  first  two  of  these  suggestions 
would  be  relatively  straightforward  (but  only  if  each  wire  remains 
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parallel  to  one  of  the  coordinate  axes)  and  would  primarily  require 
extensive  logic  for  their  numerical  implementation. 

It  should  be  noted,  however,  that  the  consideration  of  a wire 
which  is  not  parallel  to  one  of  the  coordinate  axes  would  require 
considerable  effort.  First  of  all,  integral  equation  (2.12)  would 
have  to  be  used,  and  thus  all  none  components  of  the  dyads  would  be 
needed.  In  addition,  the  benefits  of  piecewise-sinusoidal  testing 
(namely,  the  removal  of  the  differential  operator)  are  no  longer  appli- 
cable. Probably  the  most  important  problem  this  extension  would  cause, 
is  that  the  methods  outlined  in  Chapter  4 for  evaluating  the  integral 
of  the  kernel  c^uld  no  longer  be  applied.  Careful  examination  reveals 
that  this  technique  is  crucially  dependent  upon  the  fact  that  the  wire 
is  z-dlrected. 

Finally,  consider  the  application  of  this  work  to  the  case  of 
transient  excitation  such  as  an  EMP  and  characterize  the  response  of  the 
aperture  and  its  coupling  to  the  wire  within  the  cavity,  based  on  the 
obvious  application  of  the  singularity  expansion  method  (Baum  1971) , 
which  allows  one  to  express  the  total  response  in  terms  of  few  parameters 
via,  the  natural  frequencies,  natural  modes  and  coupling  coefficients. 

Or  one  can  follow  the  direct  procedure  to  obtain  the  time  dependent 
solution  by  inverse  Fourier  transformation.  Note  that  this  would  involve 
the  use  of  the  time-harmonic  solution  at  many  angular  frequencies  w 
over  the  spectrum  of  the  excitation.  It  should  be  noted  that  for  many 
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transient  excitations  (including  EMP)  the  low  frequency  portion  of  the 
I spectrum  is  dominant.  Thus  for  such  excitations,  a quasi-static  solution 

for  the  problem  is  needed  for  small  o)  in  order  to  perform  the  inverse 
Fourier  transform  numerically. 
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